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We analytically show that mesonic bound states of confined monopoles appear inside a non- 
Abelian vortex-string in massless three-flavor QCD at large quark chemical potential fi. The ori- 
entational modes CP 2 in the internal space of a vortex is described by the low-energy effective 
world-sheet theory. Mesons of confined monopoles are dynamically generated as bound states of 
kinks by the quantum effects in the effective theory. The mass of monopoles is shown to be an 
exponentially soft scale M ~ Aexp[— c(/i/A) 2 ], with the color superconducting gap A and some 
constant c. A possible quark-monopole duality between the hadron phase and the color supercon- 
ducting phase is also discussed. 
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I. INTRODUCTION 

Understanding the confinement of quarks and gluons 
is one of the most important questions in quantum chro- 
modynamics (QCD). Although there have been several 
proposals to explain the origin of the confinement, still 
no consensus has been reached. Among others, one plau- 
sible scenario is the dual superconducting picture of the 
QCD vacuum [l|: assuming the condensation of putative 
magnetic monopoles in the QCD vacuum, the color elec- 
tric flux is squeezed between a quark and an antiquark 
so that the quark- antiquark pair is confined as a meson. 
This is similar to the situation where the magnetic flux is 
squeezed into a string in the metallic superconductor due 
to the Meissner effect. Although this scenario succeeds 
in accounting for a number of properties in the QCD 
vacuum (see, e.g., [2, 3]) and is shown to be realized in 
the Af = 2 supersymmetric (SUSY) QCD 0], the con- 
densation or even the existence of monopoles cannot be 
justified in real QCD without dramatic assumptions [f|. 

If monopoles indeed exist within the theory of QCD, it 
is natural to expect that monopoles would also show up 
in QCD at finite temperature T and finite quark chem- 
ical potential fi. In the quark-gluon plasma phase at 
high T, several instances of evidence of the existence of 
monopoles and their important roles are suggested in the 
model calculations in conjunction with the lattice QCD 
simulations (for reviews, see 0,0)- One can question 
the existence of monopoles in QCD at large /i, as first 
addressed by the present authors Q. It is indeed an 
ideal situation to investigate this question at asymptotic 
large /i; the ground state is found to be the most symmet- 
ric three-flavor color superconductivity called the color- 
flavor locked (CFL) phase [9j due to the condensation 
of quark-quark pairing (for a recent review, see [13]); the 
physics is under theoretical control in this regime because 
the QCD coupling constant g s is weak according to the 
asymptotic freedom. 

In the CFL phase, the U(1)b symmetry is sponta- 



neously broken by the condensation of quark-quark pair- 
ing. This gives rise to the emergence of Abelian vortices 
(superfluid vortices) characterized by the first homotopy 
group 7ti[U(1)b] = Z [Til] - Moreover, owing to the color- 
flavor locking structure of the pairing, there als o ap pear 
non- Abelian vortices (semi-superfluid vortices) [12[ hav- 
ing only winding number 1/3 inside U(1)b and carry- 
ing a color magnetic flux. The non- Abelian vortices are 
defined as those characterized by the homotopy group 
-KiiG/H) = % for the symmetry breaking pattern G — > H 
with the condition that H is non-Abelian. The distinct 
property of non-Abelian vortices is that they have inter- 
nal collective coordinates (called the orientational modes 
or the moduli) as a consequence of the symmetry break- 
ing in the presence of each vortex. In the case of the 
CFL phase, the moduli of a non-Abelian vortex is the 
projective complex space CP 2 ~ SU(3)/[SU(2) x U(l)] 
[13]. Based on the philosophy of the effective theory (see, 
e-g-, [3 for a review) the low-energy effective world-sheet 
theory for these orientational modes near the critical tem- 
perature T c of the CFL phase is constructed in [ll| [lj| . 
The interaction between the CP 2 modes in the vortex 
world-sheet and gluons in the bulk has also been deter- 
mined [17] . 

Actually, such non-Abelian vortices originating from 
the color-flavor locking were first found in the N = 2 
SUSY XJ(N) QCD 0. Remarkably, in the Higgs phase 
of the Af = 2 SUSY QCD, the squark mass leads to the 
dynamical symmetry breaking pattern XJ(N) — > \J(1) N 
and supports the existence of monopoles characterized 
by t: 2 [V(N)/V(1) n } = Z N [lp-pH. In real QCD, on the 
other hand, it is shown in [8] that the strange quark 
mass m s together with the charge neutrality and /3- 
equilibrium conditions (required in the realistic dense 
matter) exhibits just the explicit symmetry breaking pat- 
tern SU(3) — > U(l) 2 and does not support the existence 
of monopoles dynamically as it should not. 

Still there is another mechanism supporting monopoles 
in real QCD at large fi mentioned in [8j in analogy with 
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the SUSY QCD, that is, the possible nonperturbative 
quantum fluctuations of the orientational modes. Such 
quantum effects are shown to generate a single confined 
monopole attached to non-Abelian vortices in the SUSY 
QCD [19j and a monopole-antimonopole meson in non- 
supersymmetric models motivated by the SUSY plM^ . 
If this is also the case in real QCD, monopoles must be 
confined due to the color Meissner effect of the color su- 
perconductivity in the Higgs phase. 

In this paper, we analytically show that mesonic bound 
states of confined monopoles appear as bound states of 
kinks in the effective world-sheet theory on non-Abelian 
vortices in massless three-flavor QCD at large /j,. The 
main difference from our previous analysis in Q is that 
here we ignore the effect of the strange quark mass m s , 
but take into account the quantum effects of the orienta- 
tional modes. In particular, we derive an exponentially 
soft mass scale of confined monopoles near T c : 



M 



= -c(M/A) 2 



(1) 



where A is the superconducting gap and c is some con- 
stant. 

The existence of mesonic bound states of confined 
monopoles in the CFL phase naturally realizes the "dual" 
of the putative dual superconducting scenario for the 
quark confinement in the hadron phase. We also point 
out the resemblance of the color-octet mesons formed by 
monopole-antimonopole pairs in the CFL phase to the 
flavor-octet mesons formed by quark-antiquark pairs in 
the hadron phase. This leads us to speculate on the idea 
of the "quark- monopole duality," i.e., the roles played 
by quarks and monopoles are interchanged between the 
hadron phase and the CFL phase. This duality, if real- 
ized, implies the condensation of monopoles in the hadron 
phase corresponding to the condensation of quarks in the 
CFL phase, and thus, embodies the dual superconduct- 
ing picture in the hadron phase. 

The paper is organized as follows. In Sec. HH we re- 
view the time-dependent Ginzburg-Landau Lagrangian 
(TDGL). In Sec. IIII1 we summarize the solution of a 
non-Abelian vortex and the construction of the effec- 
tive world-sheet theory on a non-Abelian vortex. In 
Sec. IIV[ we show that a mesonic bound state of confined 
monopoles appear on a vortex, and discuss its possible 
implications. Section[V]is devoted to conclusion and out- 
look. 



II. TIME-DEPENDENT GINZBURG-LANDAU 
LAGRANGIAN 



In this section, we review the TDGL Lagrangian at suf- 
ficiently large [i. We consider massless three-flavor QCD. 
This situation is different from the one considered in Q 
where the effect of the strange quark mass m s together 
with the charge neutrality and /3-equilibrium are taken 
into account. One may not ignore these effects in the 



realistic situation, e.g., inside the neutron stars. We will 
comment on this issue at the end of Sec. IIVBI 

Let us first introduce the order parameters of the color 
superconductivity, the diquark condensates <5>l,r- The 
diquark condensates are induced by the attractive one- 
gluon exchange and the instanton-induced interactions 
in the color antisymmetric channel according to the BCS 
mechanism [10] . In Dirac space, the Lorentz scalar (spin- 
parity + channel) is the most favorable, since it allows 
all the quarks near the Fermi surface to participate in the 
pairing coherently. T he p ositive parity state is favored by 
the instanton effects [2J |. The remaining quantum num- 
ber, the flavor, must be antisymmetrized for the pairing 
to follow the Pauli principle. Therefore, the diquark con- 
densate takes the form 



~ e abc e ijk ((q L )iC(q L ) k c ) 
~ tabce l3 k({qR)lC{q R ) k c ) 



(2) 



where i, j, k (a, b, c) are flavor (color) indices and C is the 
charge conjugation operator. The positive parity ground 
state is expressed by 



$1 



(3) 



We then construct the TDGL Lagrangian based on the 
QCD symmetry under 



G = 



SU(3) C x SU(3) F x U(l)j 



(4) 



where F can be either L or R and redundancy of the 
discrete groups are removed. Under the symmetry G, 
<&F transform as 



$ F e 2ie U c ^ F U F , 



(5) 



where e iB £ U(1) B , U c € SU(3) C , and U F € SU(3) F . 

Because of the absence of the Lorentz invariance in 
the medium, the Lagrangian respects the SO(3) spatial 
rotation. Near the critical temperature T c of the color 
superconductivity, the order parameters <&l,r are suffi- 
ciently small so that higher order terms in <&l,r are neg- 
ligible. Also, as long as we consider the long- wavelength 
and low- frequency deviation from the equilibrium, we can 
perform the derivative expansion. Up to the second or- 
der in time and space derivati ves, the TDGL Lagrangian 
invariant under G is given by [2a . [26| : 1 

C = Tr (K Vo&V <S> - K 3 V^V^) 

+ 2 01 4A « ' 
V = oTr + ft [Tr($ t $)] 2 + /3 2 Tr [(& <S>) 2 ] ,(6) 



To be precise, terms including the first derivative, e.g., 
Tr ($tX>o ( S > ) , related to the dissipation, are not forbidden by the 
QCD symmetry. Here we ignore them because they turn out to 
be irrelevant to the dynamics of a non-Abelian vortex finally [id . 
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where 2? M $ 



e and A arc the dielectric constant and the 



d»A v - d v A 

*ffs [Afii A v 

magnetic permeability, respectively, both of which we set 
unity in our previous works [s| Il5l Il6j j . 

The leading-order values of the Ginzburg-Landau (GL) 
coefficients a, /3i,2, and -Ko,3 are obtained from the weak- 
coupling calculations at large \x [25|, [26[ : 



a = 4N(p)lo, 



ft 
K 3 



P2 = 



T 

7C(3) 
7C(3) 



3^° 12U.T..I-- 2 



iV( M ) = 



(7) 



where iV(/i) = /i 2 /(27r 2 ) is the density of state at the 
Fermi surface and T c = 2 1/3 e 7 A/ 7r is the critical tem- 
perature of the CFL phase 27]. These GL coefficients 
can be derived by generalizing the computations known 
in nonrelativistic systems [28j. 

Using the GL potential Vgl with the GL coefficients 
in Eq. (|8]), one finds the most stable ground state, 



$ = diag(A, A, A), 



(8) 



where A = \J — a/ (8/3). This form of the ground state 
entangles the color and flavor rotations and is called the 
color-flavor locked phase. In the CFL phase, the symme- 
try G is broken down to 



H 



SU(3) C+F 
(Zs) 



C+F 



and the order parameter manifold is 



G/H 



SU(3) c _iT x U(l)i 



(Za) 



C-F+B 



= U(3). 



(9) 



(10) 



This is parametrized by would-be Nambu-Goldstone 
(NG) modes associated with the dynamical symmetry 
breaking of SU(3), which are eaten by the eight gluons 
by the Higgs mechanism, and a massless NG mode (re- 
ferred to as the H boson) associated with the symmetry 
breaking of U(l)s. 

By expanding <!> from the ground state ([5]), 

$ = Al 3 + yi ^ 1 3 + Vs ^ T g , (II) 



V2 



V2 



where (f>i and 0g (ip and £ a ) are real (imaginary) parts 
of fluctuations, mass spectra are obtained as 



m 2 G = 2Xg 2 A 2 K 3 



2a 
K~3 : 



4/3A 2 



(12) 



Here mc is the mass of the gluons which absorb C° by 
the Higgs mechanism, and mi and ms are the masses of 
<pi and 4>s ha the 3 ® 3* = 1 © 8 representation under the 
unbroken SU(3)c+f symmetry, respectively. 



From Eqs. © and (fl2]l. we have 

m G ~ 5sM, mi = 2to§ ~ 2A. 



(13) 



Because g s /j, ^> A at large /i, the CFL phase is a type- 
I superconductor as indicated by the Ginzburg-Landau 
parameters [29| : 



Kl,8 = < 1, 



(14) 



where we define two GL parameters corresponding to two 
coherence lengths l/m\$. Note that non-Abelian vor- 
tices can appear even in this type-I system, since their 
interactions are repulsive at large distances due to the 
exchange of the H boson [l3j]. This is in contrast to the 
case of the metallic (Abelian) superconductor where vor- 
tices can appear only in a type-II system with n > 1 
(under a suitable normalization). Non-Abelian vortices 
are the superfluid vortices and are created under a rapid 
rotation. 



III. NON-ABELIAN VORTICES 

In this section, we consider the properties of a non- 
Abelian vortex in the CFL phase for later discussions in 
Sec. lIVl The results are already obtained in our previous 
papers [1, [TBI, [Hj]. We also correct some of equations 
given in these references. 



A. Non-Abelian vortex solutions 

We first consider a non-Abelian vortex solution in the 
CFL phase. We make the standard ansatz for a static 
vortex-string configuration parallel to the x 3 direction 
(perpendicular to the x\-xi plane): 



$(r,0) - Adiag(e* e /M, g(r), g(r)) , (15) 

2 1 1\ 
3' 3' 3) ' 

(16) 



A t (r,8) = ±^[l-h(r)]d 
9 s ' 



lag 



with i, j = 1,2. This ansatz can be rewritten as 



$ = Ai 



F 



V3 



A; 



9s 

with profile functions 

F = f + 2g, G = f-g, 
and the U(3) generators 



gGT 8 , (17) 



(18) 



(19) 



T = -Ldiag(l,l,l), r 8 = -^diag(-2,l,l). (20) 
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Equations of motion for the profile functions f,g, and 
h are of the form 



(2k+ l) 2 TO? , , TOo 1 

A - L ^^-^(^- 3 )-^J/ = ' (21) 

(ft — l) 2 TO? . , TOr 1 , - 

A -H^--^ (A - 3) + ^T 9 ' (22) 

- 7 - ^ - 1) + / 2 (2^ + 1)) = 0, (23) 



with the Laplacian A = 9 2 



L <9 r , A = / 2 + 2<? 2 and 



P = ,f 2 — g 2 . These differential equations should be 
solved with the boundary conditions 



{f,g,h)-> (1,1,0) 
(f,g',h)^ (0,0,1) 



as r — > oo, 
as r — > 0. 



(24) 



Numerical solutions can be found in [T^, [T(| . 



B. Effective world-sheet theory on a vortex 

In this subsection we construct the effective world- 
sheet theory on a single non-Abelian vortex, which de- 
scribes fluctuations of the CP 2 orientational modes </>. 
We place the vortex-string along the x 3 axis, so we con- 
struct the effective action in the (t,x 3 ) coordinates by 
integrating over the X1-X2 plane. 

First of all, we take a singular gauge in which the single 
vortex configuration is expressed as 



= _l£^ (r) /| r 
g s V 3 



(25) 



(26) 



Then the general solution can be reproduced by acting 
the color-flavor locked symmetry on them: 

$(£/) -> U®*!!- 1 , Ai{U) -^UA*U-\ (27) 

where U € SU(3)c+f- This action changes only T$ with 
T ~ I3 unchanged. We define coordinates on CP 2 by 




T 8 IJ- L = 



(28) 



where <fi is a complex 3-column vector, and (O) denotes 
the traceless part of a square matrix O. The SU(3)c+f 
symmetry acts on <p from the left hand side as <j> — > U(f>. 
Taking trace of this gives a constraint 



1. 



(29) 



Since the phase of <j) is redundant, we find that repre- 
sents the homogeneous coordinates of the complex pro- 
jective space CP 2 . 



Physically, these degrees of freedom (called the orien- 
tational modes or the moduli) arise associated with the 
symmetry breaking from the SU(3)c+_f symmetry pre- 
served by the diquark condensates to [U(l) x SU(2)]c+f 
in the presence of each vortex 



SU(3) 



C+F 



[U(l) x SU(2)] 



CP 



C+F 



(30) 



The NG modes <f> € CP 2 propagate along the non- 
Abelian vortex-string. The form of the Lagrangian is 
determined solely by the SU(3)/[U(1) x SU(2)] symme- 
try, and is described by the CP 2 nonlinear sigma model. 
The effective Lagrangian consists of two parts: 



CeB — A 



(0) 



C 



(3) 



(0) 
eff 



(3) 
off 



off ~ '-eff 

efeida^Tr 
dxidx 2 TT 



-F Ql F m + K V &V°<S> 



2A 



F 3i F 3i + K 3 V 3 &V 3 <I> 



(31) 
(32) 

,(33) 



where the CP 2 orientational modes <j> are promoted to 
fields cj>(t,x 3 ) on the vortex world sheet. For gauge field 
we use the ansatz of Gorsky-Shifman-Yung [221 ]: 

A a (4>(x a )) = !^M[(^t) >9q ^t)] > a = ,3, (34) 
9s 

where functions p a {r) (a = 0,3) are undetermined at 
this stage, and will be determined below. Note that we 
need two independent functions in this ansatz due to the 
absence of the Lorentz invariance in medium unlike [22j , 
which was not considered in our previous paper (l6| . 
For later use, it is convenient to define the function T a 

by m 



F a {a, b) = a<j)d a <ft + bd a <$>$ + (a — b)<j>$ d a cj>$ , (35) 
with a, b € C, which satisfies 

Tr [T a (a,b)^ a (a,b)] = (\a\ 2 + \b\ 2 )£ CP 2. (36) 

Here Ccp 2 is the form of the CP 2 nonlinear sigma model 
Lagrangian: 



CP 2 



= ^d a <t> + (tfdcMtf &*</>). (37) 



By using the function T ai parts of the Lagrangian can 
be rewritten as 

v a $ = Ae* e / 3 T a (f-g + p a g,f-g-p a f), (38) 
1 



Frvi — 



5> 



^ r 2 



h(l - p a )Fa(l, 1) 



3 J 



gs 



p' a F a (l,-l). 

(39) 



By using these expressions, we calculate each term in the 
Lagrangian. First, the term of the gauge field strength 
can be calculated as 



5 



TV [F ia F ia ] = 



' ' C ^ 2 ^ - ' ^ 2ie Jl^ h{1 - Pa)p > a g a (l,i|i,-i) 



4 -/i 2 5 (l-p a )^(l,l|l,l)- , 



1 a;^ 1 

2 2~~ 



(p'a 



C 



(a) 
CP 2 ' 



( P ' a ) 2 g a (i,-i|i,-i) 



(40) 



with Q a {k, l\m, n) = Trk7-" a (/c, l)!F a {m, n)] and no sum- 
mation is taken for a. Here we have used the following 
relations 

g a (i, i|i, i) = -g a (i, -i|i, -i) = 24p 2 - (4i) 

Similarly the term including $ can be calculated to give 



Tr [D a &V a $] = 2A 2 



(1 -Pa)(f ~9? 



f(f Z +9 2 ) 



A<x) 
'-CP 2 ' 



(42) 



Substituting Eqs. (gDJ) and (g2]) into dSJ), we finally 
obtain the CP 2 La grangian 



(43) 



with two different coefficients for time and space compo- 
nents, 



SX [(p0Y + -n(l-P0Y 



2 



K 3 
An 
Ml 



(i-po)(/-.g) 2 + 



f + 9 2 



-Po 



dr- 



(p' 3 ) 2 + -(l-p 3 f 



+m 2 G [(l- p3 )(f-gf + 



f 2 +9 2 



-Pa 



,(44) 



(45) 



Co and C3 should be determined by minimizing them 
through £0,3: 

Po + y + (l-Po)^ 

-^mU(/ 2 +.9 2 )P0-(/-.g) 2 ]=0, (46) 



ft 2 

,2 



p; 3 ' + ^ + (i-p 3 )- 

-m^[(/ 2 + .g> 3 -(/-. 9 ) 2 ] =0. (47) 
From Eqs. (H3]), HU), and <gS]), C ,3 are estimated as 



C ,3 ~ 



(48) 



2 The coefficients CKo,3 in Eq. (3.18) of Ref. [TJ] should be cor- 
rected as Co, 3- 



The velocity of the CP 2 modes propagating along the 
vortex-string is then 



C3/C0. 



(49) 



This nontrivially depends on p, which we do not discuss 
in detail in this paper. By rescaling 



t i/C t, x 6 ->■ VC 3 x 6 , (50) 

the Lagrangian (|43|) is cast in the Lorentz invariant form 

£ off = £ cp2 = d a ^d a cj> + {tfd a <t>)($d a <p). (51) 

Note that the first derivative terms ignored in Eq. © 
do not contribute to the effective theory because of the 
property tr(J" a ) = [Tf|. 



IV. CONFINED MONOPOLES 

In this section, we show that mesonic bound states 
of confined monopoles appear inside the non-Abelian 
vortices by solving the effective world-sheet theory con- 
structed in the previous section in the large- N c limit. We 
further argue a possible "quark-monopole duality" be- 
tween the hadron phase and the color superconducting 
phase. 



A. Bound state of monopole-antimonopole: 
Kink-antikink pairing on a vortex 

In this subsection, we consider the properties of the 
solution to the CP 2 nonlinear sigma model Lagrangian 
for the orientational modes, by taking into account the 
quantum effects. Thereby we will find that there appear 
a kink-antikink pairing on a vortex which can be iden- 
tified as the mesonic bound state of a monopole and an 
antimonopole in the 3+1 dimensions. 

Unfortunately, however, the solution to the CP 2 non- 
linear sigma model is not known so far, although the CP 1 
model [equivalent to the 0(3) nonlinear sigma model] 
is solved rigorously 30]. Here we consider the CP^ -1 
model instead and solve the model to leading order of 
1 /N c following [3l], [HJ . Owing to the qualitative similar- 
ity of the solutions to the CP 1 and CP Na ~ 1 models, the 
solution to the CP 2 model should be approximately de- 
scribed by the solution to the CP Wc_1 model with taking 
N c = 3 at the end. This is the only assumption which 
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we will make in our calculations. To make our paper 
self-contained, we shall describe the original derivation 
[3ll I32I ] for our Lagrangian (fSTj) under the constraint (|29|) 
in the following. 

We first perform the Hubbard-Stratonovich transfor- 
mation by introducing the auxiliary field to eliminate the 
quartic term in Eq. (|51[) . Because the quartic term is the 
vector-vector type interaction, the auxiliary field should 
be the gauge field A a . After the Hubbard-Stratonovich 
transformation, the Lagrangian becomes 



C eff = (d a - iA a )tf (d a + iA a 



(52) 



with the constraint (|29[) . Actually, eliminating A a by us- 
ing the equation of motion for A a , one can easily check 
that Eq. (|5"2"|) reduces to Eq. (|5"Tj) . Equation (|5"2"j) can 
be regarded as the U(1)d gauge theory ("D" stands 
for a dummy gauge symmetry); it has the local gauge 
symmetry <fi — » e ^ x '4> under the gauge transformation 
A a — > Aa — d a 8(x). To take into account the constraint 
(I29p in the Lagrangian, we then introduce another aux- 
iliary field (7 as a Lagrange multiplier, to obtain 

£ c ff = {d a - iA a )<ft (3 a + iA a )^> -a^cp-l). (53) 

This expression of the CP^^ 1 model is nothing but the 
Kahler quotient. After rescaling the </> and (jr variables, 



CoC*3 



1/4 



C0C3 



1/4 



4>\ (54) 



the partition function of the theory is given by 



UadA a ] e tS \ 



iA° 



Z = 

where the action is given by 
S = I dxodxs 



(55) 



(d a -iA a )^(d a 
N, 



V CqCj 



(56) 



Note here that the coefficient of the kinetic term becomes 
unity due to the rescaling of the fields, Eq. (j54"|) . together 
with the rescaling of measure, Eq. (l50|) . The prefactor 
of N c /3 in Eq. (|S"6"|) is introduced to perform the 1/N C 
expansion consistently below, and is chosen such that it 
reduces to unity for N c = 3. 

Integrating out (j> and (ft , one obtains 



Z = J [dadA a ] exp 



iV c trln -{d a + iA a Y - a 



N r 

+i^\/CoQi / dx dx 3 a 



(57) 



The Lorentz invariance implies the saddle point A a = 
and constant a. 



To leading order of 1/N C , varying the partition func- 
tion with respect to a gives the gap equation: 



. v 7 C0C3 



A=A 



d 2 k 



1 



(2tt) 2 k 2 



0, 



(58) 



where we have introduced the cutoff of the low-energy 
effective theory, A = A. This is not a dynamical cutoff 
of the CP 2 model but is a physical cutoff, namely A is the 
mass gap of the quasiparticles of quarks in the original 
GL Lagrangian [see Eq. (fl3l) ]. After the integral, one 
arrives at 



M 



Using (j48|) . M is expressed as 

M - Ae" c(Al/A) 



(59) 



(60) 



with some constant c. M can be now identified as the 
mass of cf> and <fr induced by the quantum effects from 
Eq. (|5"6"|) . This mass gap for the orientational modes is 
required by the Coleman- Mcrmin- Wagner theorem in the 
1+1 dimensions, as mentioned in [8J. 

We then consider the fluctuations around the saddle 
point A a = and a given in Eq. (|59|) . For this pur- 
pose, we expand the partition function with respect to a 
and A a . Since higher order terms in a and A a are sup- 
pressed in the large- N c limit, only the quadratic terms 
are relevant in the following. Expansion of the functional 
determinant in Eq. (|57p can be understood in terms of 
the Feynman diagrams. It turns out that the relevant 
diagrams are the propagator of A a at one-loop level 
[the U(1)d "photon" self-energy]. Finally, the dynami- 
cally generated kinetic term of the gauge field reduces to 
cNa{-g^ v k 2 + k^kv) with the coefficient c = 1/(12ttM 2 ) 



[32|. The form of the kinetic term is fixed by the gauge 
invariance. 

Now the effective world-sheet theory including the 
quantum effects to leading order of 1/A C is summarized 
as 



>C^ ant = (d a -iA a )tf(d a + iA a ) 

N c F 2 

~48ttM 2 q/3 ' 



M' 



(61) 



By rescaling A a so that the kinetic term of A a is canon- 
ically normalized, 



\2ttM 2 a 

A a ^J^-A 



the effective Lagrangian reduces to 



(62) 



CT nt = [da- iAaMj 1 -^-) (V + iA a RI^J-^— 



1 



-My - -F 2 a[j . 



(63) 
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This implies that (f) an d ft have the effective U(1).d 
charges ±My / 12n/N c . Since we are considering the 1+1 
dimensions, (j) & n d (ft 1 * are confined by the linear potential 



V{x,y) = — — |a;-y| 



(64) 



where x and y are the a; 3 coordinates of (j> and (ft . 

We are now ready to understand the confining poten- 
tial between <j> and (ft from the 3+1 dimensional view- 
point. Remembering that <f> is the orientational moduli 
of the non-Abelian vortex, a quantum state of CP 0-1 
model is in one-to-one correspondence with a quantum 
vortex state. Since there exists only one ground state in 
the CP^ -1 model, so is the quantum vortex state whose 
orientation is not fixed in a particular direction. When 
4> and (ft are placed at positions x and y, because of the 
linear potential between them, the string tension of the 
vortex between x and y is larger than that outside this 
region by ~ M 2 /N c ; the inner vortex is an excited state 
compared with the outer vortices (ground state) . 

There is another perspective for understanding of this 
phenomenon (22|. The vacuum structure of the CP^ -1 
model can be realized by looking at the ^-dependence of 
the theory [H [H Q 



(65) 



where the gauge field A a is the one after rescaling (|62|) . 
Recalling that the vacuum energy E(9) is of order N c in 
the large-iV c limit, E(9) is expressed as 



E{9) = NJ 



(66) 



Here f(9) is an even function of 9 due to the CP sym- 
metry under which 9 transforms as 9 — > —9. E{9) must 
also satisfy the periodicity, 



E{9) = E(9 + 2ir). 



(67) 



One might suspect that these two conditions arc incom- 
patible at first sight. However, there is a way out: both 
of them can be satisfied when E{9) is a multibranched 
function as 



E{9) = N c min / 

k 



9 + 2nk 



(fc = 0,l,.- - ,N C -1). 



(68) 



Expanding f{9) = /(0) + (l/2)/"(0)6> 2 + • • • and consid- 
ering that higher order terms in 9 are suppressed at large 
N c in Eq. (|68| . the vacuum energy at 9 = is given by 



M 2 

E(0)=E + C—k 2 , (fc = 0,l r 



.Nr. 



1 



(69) 



with some constant C. Therefore, there exist N c local 
minima among which only one is a true ground state 
while the others are quasivacua, see Fig. Q] 



It is now natural to interpret <f> and (ft as a kink and 
an antikink interpolating the adjacent local minima on 
a vortex, respectively [221 ] ; taking into account the codi- 
mension, this bound state neutral to the XJ(1)d charge 
can be identified as the bound state of a monopole and 
an antimonopole in terms of the original 3+1 dimensions, 
as illustrated in Fig. [TJ a monopole and an antimonopole 
with the mass M are confined into the mesonic bound 
state by the linear potential. A similar understanding 
has been demonstrated in [35[ based on the comparison 
with the SUSY QCD. 





FIG. 1: A schematic illustration of the nonperturbative po- 
tential and kinks interpolating between the ground state and 
the metastable states, in the cases of (a) N c = 2 and (b) 
iV c = 3. Kinks {<p and ft) can be identified with monopoles 
(M and M) from the bulk 3 + 1 dimensional point of view. 
The total configuration is a bound state of a monopole and 
an antimonopole. 



Let us discuss the representations of these objects. 
The fields cf> and (ft in the effective theory transform 
as N c and N c * (anti)fundamental representations un- 
der SU(A f c ), respectively. The degrees of freedom of <p 
is 2{N C — 1) after fixing U(1)d gauge symmetry. For 

instance 6 = , 1 (l,6i) T for N c — 2, and 6 — 
\A+I h il 2 

1 =(1, b 1 ,b 2 ) T for N c - 3. For A^ c = 2, (0t) 



y/l + \bi\ 2 + \b2\ 

represents one (anti)kink, as can be seen in Fig. [TJa). 
Each of them corresponds to one (anti)monopole in the 
bulk. For A^ c = 3, which is the case of the CFL phase, 
one (anti)monopole is a composite state of N c — 1 = 2 
(anti)kinks, each of which has one complex moduli (po- 
sition and phase), as seen in Fig. [rfb). 
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We conclude that (anti)monopoles belong to N c (N c *) 
fundamental representations of SU(iV c ), and they appear 
as a mesonic bound state. This mesonic bound state 
belong to N c ® N c * = 10 — 1 representation. It was 
shown in [30( that the singlet in this decomposition does 
not appear in the spectrum in the CP 1 model (N c = 2). 
This was interpreted in [35j that the singlet corresponds 
to a set of monopole and antimonopole with opposite 
charges, which is unstable to decay. Although there is no 
such a calculation for N c > 3, we expect that the same 
holds. 

Before closing this subsection, let us make one com- 
ment on fermions. As discussed in (32l. [3t|. fermions can 
be incorporated to the CP^" 1 model. In fact, quasipar- 
ticles of quarks are shown to be trapped in the core of 
non-Abelian vortices in [371 ] , in which fermion zero modes 
belonging to the triplet of the SU(2) unbroken symme- 
try in the core of the vortex have been found. However 
coupling to the bosonic CP 2 model is not known yet. 

In summary, we found a mesonic bound state of con- 
fined monopoles with the mass M given in Eq. (1601) in- 
duced by the quantum effects inside a non-Abelian vor- 
tex. 



B. A possible quark- monopole duality 

In this subsection, we would like to ask the implica- 
tions of our results. In Sec. IIV A[ we found the color- 
octet mesonic bound states formed by of monopolc- 
antimonopole pairs. Because of the color-flavor lock- 
ing, they also form flavor-octet under the remaining 
SU(3)c+f- Clearly, these bound states resemble the 
flavor-octet mesons formed by quark- antiquark pairs in 
the hadron phase. This leads us to speculate on the 
idea of the "quark-monopole duality": the roles played 
by quarks and monopoles are interchanged between the 
hadron phase and the CFL phase. If this is indeed the 
case, this would imply the condensation of monopoles in 
the hadron phase corresponding to the condensation of 
quarks in the CFL phase. This naturally embodies the 
dual superconducting scenario for the quark confinement 
in the hadron phase [1| . 

The possible quark-monopole duality may have some 
relevance to the one-to-one correspondence of the physics 
without any phase transition between the hadron phase 
and the CFL phase conjectured by Schafer and Wilczek 
38]. This is called the "hadron-quark continuity" and 
may be realized in the QCD phase structure in three- 
flavor limit as explicitly shown in [39|, H3|. 3 The cor- 
respondence in the quark-monopole duality and the 
hadron-quark continuity is summarized in Table [I] The 



idea of the hadron-quark continuity is supported by a 
number of nontrivial evidences: the same symmetry 
breaking patterns, the fact that confinement phase is in- 
distinguishable from the Higgs phase 41], the one-to-one 
correspondence of the elementary excitations such as the 
baryons, vector mesons [42j], and pions [1(3], 4 and the 
equivalence of the form of the partition functions in a fi- 
nite volume called the e-regime [43| , between the hadron 
phase and the CFL phase. 



The quark-monopole duality raises a question regard- 
ing possible other states formed by monopoles. In the 
hadron (confining) phase, a baryonic bound state is made 
of three quarks. It has been found in [44j by the lattice 
QCD simulations that three quarks are connected by a 
Y-shaped junction of color electric flux tubes. What is 
the counterpart in the CFL phase? We expect that it is 
a junction of three non-Abelian vortices with total color 
fluxes canceled out at the junction point: red, blue, and 
green color magnetic flux tubes meet at one point, see 
Fig. [2j We note that they carry correct the baryon num- 
ber as we expect for a baryon; each flux tube carries the 
U(l) b winding number 1 /3, and all of them join together 
to constitute one U(1)b vortex with the U(1)b winding 
number one. However we have not specified the electro- 
magnetic charges of fluxes at this stage because we have 
ignored the electromagnetic coupling of vortices. A simi- 
lar string junction (without monopoles) is known to exist 
in a U(l) x U(l) model [45| . The configuration in Fig. [2] 
cannot be discussed in the effective field theory of a sin- 
gle vortex anymore, but one may be able to do that by 
considering multivortex effective theory. 5 



3 Here we mean the "hadron phase" by the three-flavor symmetric 
nuclear matter (hyper nuclear matter) where the U(1)b symme- 
try is dynamically broken by the baryon-baryon pairing. 



4 The difference of the singlet can be ascribed to the mass splitting 
between the octet and singlet. For the quarks in the CFL phase, 
the singlet is twice heavier than the octet [see Eq. Ql)], which 
is expected to correspond to the excited singlet baryonic state in 
the hadron phase. For the vector mesons in the hadron phase, 
the mass splitting is induced by the diquark condensate, and 
the flavor singlet vector disappears at some intermediate fi [4 J . 
For the NG modes, the singlet r/ meson in the hadron phase is 
heavy due to the U(1) J 4 anomaly, but becomes a light NG mode 
in the CFL phase by the instanton suppression at large /i |24j . 
Therefore, the hadron-quark continuity still works. 

5 Multivortex states were studied in the SUSY QCD [4rj |. in which 
case no static interactions exist between vortices when they are 
placed parallel to each other. In our case of the CFL phase, par- 
allel vortices are repulsive at least when they are well separated 
[T3 |. However short range interactions have not been studied 
yet, and there is a possibility of attraction at short distance. In 
any case, we consider that the bound state should be quantum 
mechanically (but not necessary classically) stable with the ap- 
pearance of monopoles, as a meson of monopoles found in this 
paper; three vortices with different color fluxes join to one V(1)b 
vortex with no fluxes. 
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Phases 


Hadron phase (hyper nuclear matter) 


Color-flavor locked phase 




Confinement 


Higgs 


Quarks 


Confined 


Condensed 


Monopoles 


Condensed? 


Confined 


Coupling constant 


Strong 


Weak 


Order parameters 


Chiral condensate (qq) 


Diquark condensate (qq) 


Symmetry breaking patterns 


SU(3)l x SU(3) fl x U(1) B 


SU(3)c x SU(3) L x SU(3) B x U(1)b 




-> SU(3) i+fl 


-)• SU(3) C +L+R 


Fermions 


Octet baryons 


Octet + singlet quarks 


Vectors 


Octet + singlet vector mesons 


Octet gluons 


Nambu-Goldstone modes 


Octet pions (qq) 


Octet + singlet pions (qqqq) 




H boson 


H boson 



TABLE I: Comparisons of the physics between the hadron phase and the CFL phase in massless three-flavor QCD: symmetry 
breaking patterns [the U(1)a and discrete symmetries are suppressed here] and the elementary excitations. We have shown 
that confined monopoles appear in the form of the mesonic bound state in the CFL phase. Still one missing piece in the table 
is the properties of monopoles in the hadron phase, for which we speculate the condensation of monopoles corresponding to 
the condensation of quarks in the CFL phase. See the text for further explanations. 



hadron-quark continuity survives when one turns on m s ; 
there are other candidates for the ground state at inter- 
mediate fj, other than the CFL phase under the stress of 
m s , such as the meson condensed phase, the crystalline 
Fulde-Ferrell-Larkin-Ovchinikov phase, gluon condensed 
phase, etc fioj ] . 



CONCLUSION AND OUTLOOK 



FIG. 2: Baryonic bound state of three monopoles. Three 
monopoles are connected by a Y-shaped junction of color 
magnetic flux tubes with the total color canceled out. U(1)b 
vortex is represented by <g> at the junction point. This state 
is dual to a baryon made of three quarks connected by a Y- 
shaped junction of three color electric flux tubes in the hadron 
(confining) phase. 



Finally, let us note the effect of the strange quark mass 
m s with the charge neutrality and the /3-equilibrium con- 
ditions, as is expected in the physical dense matter like 
inside the neutron stars. This situation is considered 
previously without the quantum effects for the orienta- 
tional modes and the nonexistence of monopoles in the 
CFL phase is shown Q. Even if we further take into 
account the quantum effects, they are negligibly-small: 
the scale of the potential m 2 s /g s for the orientational 
modes induced by m s [8|, is much larger than that in- 
duced by the quantum effects M 2 , for realistic values of 
the parameters, m s ~ 100 MeV, fi ~ 500 MeV, and 
A ~ 50 MeV; confined monopoles will be washed out 
by m s . Therefore, we expect that the notion of the 
quark-monopole duality is well-defined close to three- 
flavor limit. It is also a dynamical question whether the 



In this paper, we have analytically shown that mesonic 
bound states of confined monopoles appear in the color- 
flavor locked phase of three-flavor QCD at large quark 
chemical potential fi. They are dynamically generated as 
kinks by the quantum fluctuations in the effective world- 
sheet theory for the orientational modes on a non-Abelian 
vortex. The mass of monopoles has been computed as 
M <~ Aexp[— c(fi/A) 2 ] with the superconducting gap A 
and some constant c. 

Both of the mesonic and baryonic bound states of 
monopoles studied in this paper have long fluxes extend- 
ing to spatial infinity. In a realistic situation these may 
not be appropriate because they have infinite energy. In 
order for them to have finite energy, such long fluxes can 
be made as loops, as illustrated in Fig. [3] 6 For a me- 
son, one can check if this can occur by studying the CP 2 
model on compactified space S . For a baryon the situ- 
ation would be more difficult. A possible configuration 
for a baryon is given in Fig. [3Jb) . 



A similar configuration to Fig. [3j a) is also discussed in the QCD 



vacuum 



6] 
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(b) 

FIG. 3: Illustrations of (a) a meson and (b) a baryon with 
finite energy. Flux tubes make loops, (a) corresponds to a 
vortex loop with monopole-antimonopole, which looks like a 
necklace. In (b) three color magnetic fluxes joint at two ends 
of a single U(1)b vortex denoted by a broken line. 



Before closing the paper, let us address several impor- 
tant questions to be investigated in the future. 

1. Although we have shown that confined monopoles 
dynamically appear as the monopole-antimonopole 
mesons by the quantum fluctuations, their topo- 
logical properties are still unclear. First of all, 
one should clarify the homotopy group responsi- 
ble for the existence or the topological stability of 
monopoles. One should also calculate the color 
magnetic charge of the confined monopoles and 
check if the Dirac condition for the color mag- 
netic charge is indeed satisfied. The color magnetic 
flux should be related to the dynamically induced 
V(1)d gauge field. These are in contrast to the sit- 
uation in the SUSY QCD: flux matching between 
a monopole and a vortex attached to it has been 
demonstrated |47j |. 



2. Our derivation is based on the effective world- 
sheet theory on a non-Abelian vortex derived from 
the time-dependent Ginzburg-Landau Lagrangian. 
This is only valid near the critical temperature of 
the CFL phase. One should argue the existence 
of monopoles far away from the critical tempera- 
ture, e.g., at T = 0. As this concerns, the analysis 
beyond the Ginzburg-Landau theory can be stud- 
ied by the Bogoliubov-de Genne equations which 
describe condensates and quasiparticles from the 
fermion degrees of freedom. In fact fermion modes 
have been studied in the presence of a non-Abelian 
vortex by the Bogoliubov-de Genne equations [37| ■ 

3. In the case of the SUSY QCD, quantum effects in 
the (1+1)D on a vortex can be explained by instan- 
ton effects in the original (3+l)D theory [l9l. l2lj. 
Actually, instantons can stably exist inside the 
vortex world-sheet [4§ |. In real QCD at asymp- 
totic large fj,, bulk instanton effects with the en- 
ergy ~ l/g 2 3> 1 are highly suppressed due to the 
asymptotic freedom of QCD and the screening of 
instantons [24] . The instanton energy in the vortex 
world-sheet is C , 3 ~ (m/A) 2 > l/g 2 [see Eq. (l4"g|) ]. 
which is further suppressed, consistent with our re- 
sult. As in the SUSY QCD, the quantum effects 
inside the vortex may be explained by instantons 
trapped in it, which remains as a future problem. 
On the other hand, the fact that the instanton en- 
ergy Co, 3 ~ (^/A) 2 inside the vortex is larger than 
the one ~ l/<? 2 in the bulk implies that instantons 
are repulsive from vortices. 

Note added. — While this work was being completed, we 
learned that A. Gorsky, M. Shifman, and A. Yung [49| 
have independently found very closely related results. 
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